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The lightest supersymmetric particle (LSP) is stable in an i?-parity conserving 
theory. In this article the steps needed to calculate the present day mass density 
of such a particle are detailed. It is shown that there can be a significant amount 
of LSP dark matter in the universe. Furthermore, relic abundance considerations 
put an upper bound on how large supersymmetry breaking masses can be without 
resorting to finetuning arguments. 



1 Introduction 

The most general gauge invariant superpotcntial will generally lead to unac- 
ceptable fast proton decay. For this reason, a discrete symmetry must be 
posited that banishes baryon and/or lepton violating operators whick con- 
tribute to this decay. The simplest such discrete symmetry is i?-paritya. Ex- 
act i?-parity conservation makes the lightest supersymmetric partner (LSP) 
stable, thereby introducing many interesting phenomena. The most studied 
phenomenon is the large missing energy signature associated with production 
of superpartners and subsequent decay into the LSP plus jots, photons, lep- 
tons, etc. For experimental reasonsH and theoretical reasonstl the LSP is now 
expected to be the lightest neutralino. 

A stable LSP has more than collider physics consequences. If they are 
created in the hot and violent early days of the universe, then there should 
be some left over today, and perhaps they could have significant cosmologi- 
cal and astrophysical consequences. If it turns out that the LSP constitutes 
a significant mass fraction of the universe then it should be possible to wit- 
ness large-scale gravitational effects of these particles on galaxies and clusters. 
Experiments have been testing large-scale gravity for decades now, and the 
current consensus states that there are non-luminous sources of gravitational 
import beyond the ordinary baryonic matter that makes up planets and starsu. 
Part of the evidence of additional mass-energy beyond our luminous matter 
includes rotation curves of galaxies and infall of clusters. The evidence for 
non-baryonic dark matter can be attributed to the successful agreement be- 
tween measurement and theory in big bang nucleosynthesis (BBN). BBNJjcUs 
us that baryonic mass fraction of the universe is probably less than 10% (I 
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am making the usual assumption that the total energy density of the universe 
is the critical energy density in accord with a general inflation scenario.) 

One is left wondering what the rest of the universe is made of. One sug- 
gestion is a weakly interacting massive object, or WIMP. The acronym WIMP 
is a good general label, but it is a bit misleading. Upon closer inspection a 
particle of 0{mw) mass which has full-strength SU{2) (weak) interactions is 
generally not a good dark matter candidate, and yields a relic abundance less 
than is required to have astrophysical significance. Some additional suppres- 
sions are generally needed in the annihilation cross-section. As we will see in 
subsequent sections, the LSP generally has the right suppressions to make it 
a good dark matter candidate while at the same time having mass near the 
weak-scale. For this reason, I will not use the generic word WIMP, and instead 
refer to the dark matter candidate as the LSP. LSPs are excellent candidates 
for the dark matter because if enough can he around they allow conformance 
to experimentally determined properties of galactic rotation curves, structure 
formation and big bang nucleosynthesis. Furthermore, calculation of their relic 
abundance indicates that LSPs could contribute most of the energy density of 
the universe. This is a non-trivial separate test that confirms that enough can 
he around to solve the observational problems. 

In the lightest neutralino, supersymmetry provides a natural dark matter 
candidate. In fits of optimism one could even declare the galactic rotation 
curve anomalies, etc. as positive experimental evidences for supersymmetry. 
However, in this chapter the focus will be on two main topics. First, and 
foremost, I will provide the details on how to calculate the relic abundance 
of a weakly interacting particle. Since obtaining the correct relic abundance 
is a somewhat involved calculation, but also an important one, it is useful to 
have detailed discussion. Second, I will tailor some additional remarks about 
the calculation to the lightest neutralino, and show how the results constrain 
other supersymmetric particle masses. In general it can be shown that there is 
an upper limit to superpartner masses due to relic abundance considerations 
alone. The limit is based on physical necessity and not on arbitrary fine- 
tuning considerations. This insight is as important as the realization that 
supersymmetry can cure the "dark matter problem." 

The subsequent sections reflect the goals presented in the previous para- 
graph. Much of the emphasis will be placed in the techniques of calculating 
the relic abundance. I have attempted to include in this one source all the nec- 
essary general relativity, statistical mechanics, and particle physics knowledge 
needed to follow a precise calculation. I have also included a section which de- 
rives an accurate and approximate solution to the Boltzmann equation. These 
results will then be used to analyze quantitatively how the supersymmetric 
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spectrum is affected by the rehc abundance constraint. 

2 Solving the Boltzmann equation 

The starting point is the Boltzmann equation, 
(Iti 

— ^-mn-{av){n'-n%) (1) 

where H is the Hubble constant, n is the particle number density in question, 
Ueq is the particle number equilibrium density, and {av) is the thermal aver- 
aged cross-section. The Boltzmann equation is simple in form but somewhat 
subtle to solve. The differentiating parameter is time t, however n and Tigg are 
most easily characterized by temperature. Furthermore, the Hubble constant 
evolution is best traced by the relative time change in the scale parameter 
H — a/a. Only one parameter is independent, and so the first step will be 
to cast the Boltzmann equation into a purely temperature dependent relation. 
The two equations that will allow us to do that are the "Friedmann equation" 
and the "conservation of entropy equation." 

The use of Einstein's equation of general relativity is necessary to reveal 
the explicit time dependence of a/ a. This familiar equation states that 

R^lu — -^Rgfj-v = SnGNTf^i,. (2) 

To do explicit calculations the metric and stress-energy must be defined. We 
use the flat Robertson- Walker metric 

ds^ ^dt^ -a'^{t)[dx^ + dy^ + dz^] (3) 

which assumes the universe to be fiat, homogeneous, and isotropic. This is 
equivalent to the metric tensor 

g^, = diag(l, -a\t), ^a^{t), -~a^{t)). (4) 

The general stress-energy consistent with an homogeneous and isotropic uni- 
verse is 

T^^ = diag(p, -p, -p, ~p). (5) 

Einstein's equation is of course valid for each {/^t^}, however jwe only need the 
00 component. From the metric tensor it is straightforward Q to compute the 
Ricci tensor component i?oo and the Ricci scalar R: 



Roa — 



(6) 



(Dots indicate time derivative.) The 00 component of the Einstein equation 
under the flat Robertson- Walker metric is then simply, 



- = — where k = W " ■ (8) 
an V ottGn 

This equation is often called the Friedmann equation. This will be quite useful 
to get rid of the scale factor in the Boltzmann equation. The value of p on 
the right hand side of the equation will be calculated later and is conveniently 
parameterized by its temperature dependence. Thus, the Friedmann equation 
provides a nice connection between the time dependent relative scale factor 
(the Hubble constant) and the temperature. 

In thermal equilibrium entropy is conserved, and using the first law of 
thermodynamics we can identify the entropy as S{T) = {p + p)V/T up to an 
irrelevant constant. With V = we define the entropy density s{T) as 

The conservation of entropy means that S{T) = s{T)a^ is time independent: 
^(s(Ty)^O ^ t^-^l§^y (10) 

(The prime on s'{T) indicates a temperature derivative.) This equation is a 
direct result of the conservation of entropy in thermal equilibrium and so is 
called the "conservation of entropy equation" . Its utility is relating the time 
derivative of the temperature (T) to the scale factor time derivative (a/a). 

We actually have enough information from the above paragraphs to con- 
struct the temperature dependent Boltzmann equation. First, we rewrite 
dn/dt = {dn/dT)T and use the "conservation of entropy" equation to replace 
T in favor oi a/a. Then we use the "Friedmann equation" to replace d/a in 
favor of the energy density p{T). The result is 

dn s'(T) [ kJU) , o o ,^,,1 

where J(T) = {av){T) is the thermal averaged cross-section. 

This equation might not look like much progress; however, all non-trivial 
dependences from p{T) and s{T) are easily calculated as functions of tem- 
perature. This will be demonstrated below. At sufficiently high temperature 
(Tff), where all relevant particles are in thermal equilibrium, then we know as 
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a boundary condition that niTn) = ngqlTn)- One then need only integrate 

down to today's temperature (T ~ 0) to obtain the current number density in 
the universe n(0). The mass density is then just = m^n(0) where rriy^ is 
the mass of the reUc particle of interest. 

Often it is of interest to compare the mass density of our relic particle to 
the critical density pc = k^H^ needed for a flat universe. The relevant formula 
is 

= (-) 

For a flat univc!rse the; smn of all contributing fi's (baryons, neutrinos, cos- 
mological constant, ncutralinos, etc.) must be equal to 1. If a massive stable 
particle makes up a significant fraction of the total critical density then it is an 
interesting cold dark matter candidate. If the calculated critical density is too 
high, then it is said to "overdose the universe" , meaning that the universe be- 
came matter dominated too early and it is impossible to reconcile the current 
mass density and the Hubble constant with the age of the universe. 

In order to effectively solve the Boltzmann equation we must have an un- 
derstanding of the thermodynamic quantities which enter the equation. These 
are the equilibrium number density neq{T), the energy density p{T), and the 
entropy density s{T). Since s{T) = {p + p)/T we can focus on calculating 
the pressure p{T) rather than computing s(T) directly. Calculating thermo- 
dynamic quantities is standard statistical mechanics and can be found in nu- 
merous sources. Here I will merely argue the most salient points that will lead 
to workable equations. 

The density of states in a phase space volume d^xcfk is l/(27r)^. Integrat- 
ing over the volume, the phase space volume density of states is 

where k = \k\ here. We are interested to begin with in the number of particles 
per unit volume (number density) . It will be necessary to multiply the density 

of states times the mean occiipation mimbcr for a given momentum state |fc). 
The mean occupation number is expressed by the Fermi and Bose distribution 
functions, 

where = — 1, 1 for bosons and fermions respectively. Thus, the total number 
density integrated over all momentum modes is 

//•OO 
d^f^{k,T) = ^ dkef^{k,T) (15) 
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where g is the number of internal spin degrees of freedom. 

Similarly, we can calculate the energy density and pressure as 



/•CO 

P^^) = A / dkkWk^+m^Uk,T) 

Jo 



dk- 



(16) 
(17) 



Since s(T) = {p{T)+p{T))/T we are done. We should keep in mind that neg(T) 
in the Boltzmann equation only applies for the one relic particle. However, the 
p(T) and s{T) in the equation are the total energy density and entropy density 
summed over all particles. 

We always want T to be identified with the photon temperature. The ther- 
modynamic quantities technically should be calculated at the temperature of 
each particle Ti and then the contributions of each particle to the density and 
pressure should be summed. This creates a subtlety when a stable particle de- 
couples from the photon thermal bath yet still contributes to the mass density, 
pressure, and entropy of the universe. When a particle decouples its entropy 
density is separately conserved from that of the photon bath. However, at 
the decoupling temperature Td the sum of the two entropies for T < Td must 
be equal to the total entropy for T > Td- This is just entropy conservation. 
Mathematically this is expressed as s{Ta)T^ = s{Tb)T% where Tb = Td + 6 
before decoupling of particle x ^-nd T — = T^, and Ta = Td — 5 \s the 
temperature after decoupling and T — T-y ^ T-^. However when ^ — > 0^ we 
can identify Ta = T and Tb 



T^ to obtain 



T,, = T 



siTd - S) 



-,1/3 



siTd + S) 



(18) 



For stable particles which have decoupled then Ti should be substituted into 
the formulas for p(T) and s{T). 

When additional particles "distribution decouple" from the thermal bath 
{W, h, etc. whose density goes to zero as T <C m) the T in Eq. |l^ is no 
longer the current photon temperature but rather the photon temperature 
before the "distribution decoupling" . Using conservation of entropy again at 
the decoupling boundary, we can find the relation between and the current 
photon temperature (T = T^): 



lim T 

s-,o+ 



s{Td ~~ 5) 



siTd + S) 



1/3 



n 

T^<mi<Ti 



s{T^ + S) 



1/3 



(19) 
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_s[Td + S)_ 

The thermodynamic quantity integrals are straightforward to calculate 
numerically for any mass and any temperature. They also can be solved an- 
alytically in the limits that T ^ m and T <^ m. For example, if T ^ m 
then the boson energy density is gT^Tr^/SO, and if T <C m then energy density 
decreases rapidly as p '--^ rp3/2^-m,/T ^ Doing the same for fermions we can 
construct the following approximation for p(T) and s(T): 

PiT) = ^PiT)T' (21) 

siT) = -^s{T)T' (22) 

where 

p{T) ^ Yl 9^0m-m,)(^Y .9.em-m0 (1)23) 

i— bosons i— f erniions 

s{T) c Y fM0(T,^m,)(^y + 1 Y 9^0{T,~m,)f^^2A) 

i— bosons i— fermions 

Again, for most particles Ti — T; stable particles decoupled from the photon 
will have Ti 7^ T. The equilibrium number density is neq{T) ~ 1.2gT^ /ir'^ , 
0.9gT^ /tt^ (bosons, fermions) for T ^ m and decreases rapidly for T <^ m. 

3 Approximating the relic abundance 

Often it is useful to employ apgECffiimation techniques to calculate efficiently 
and reliably the relic abundanceQoB'Q. It is based primarily on dividing thermal 
history into two distinct eras: before freeze-out, and after freeze-out. Freeze- 
out means that the annihilation rate of the particle, F = neq{crv), is less than 
the Hubble expansion rate. Once this occurs the particle can no longer remain 
in equilibrium. 

From previous sections we have the tools to solve for this freeze-out tem- 
perature given the condition that F = iJ. Massive weakly interacting particles 
are non-relativistic when they freeze-out, and so the equilibrium number den- 
sity can be well approximated as 



ne,=5(^)'^%-™/^, (25) 
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where g is the number of degrees of freedom of the rehc particle (2 for a 
neutraUno) . It is convenient to reexpress all temperatures as the dimensionless 
variable x = T/m-^ such that 



9^ .2.3/2g-l/x^ 



(26) 



(27r)3/2- 

Furthermore, from our discussion earlier we know that the Hubble constant is 



H = - = 
a 



(27) 
(28) 



where A = y/STrHT/AE. The freeze-out condition that neq{(Jv) ^ H aX x ^ xp 
yields a transcendental equation for xp. 



In 



gm{(7v){xF) 



/]VFV^^(27r)3/2 



(29) 



Eq. ^ is only logarithmically dependent on (av) and yields a value of — 
1/20 quite consistently for massive weakly interacting particles. 

After freeze-out the actual number density remains high above the sub- 
sequent would-be equilibrium number density. It is therefore appropriate the 



approximation n — 
niently rewritten as 



The Boltzmann equation can then be conve- 



dx 



(30) 



AVNf 

where / — n/T^ subject to the boundary condition fixp) — n^qjTp. To find 
the current number density we must integrate this equation from xp (freeze- 
out) down to a; ~ (today): 



/(0) = 



AyfW^ 



m^Jixp) + Ay^Npf^^^xp) mxJ{xp) 



where 



J{xp) 



dx{av) (x) 



The calculated relic abundance is then simply 



J{xp) \T. 



(31) 
(32) 

(33) 
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Tp — m^xp Np 





- m-c 


494/8 


rric 


— TTIt- 


578/8 


rriT 


- mb 


606/8 


rrib - 


- mw 


690/8 


mw 


- mz 


738/8 


mz 


- mt 


762/8 


> mt 


846/8 



Table 1: Degrees of freedom corresponding to freeze-out temperature Tp = xpm^ ~ m-^/20. 



The ratio /T^ < 1 is the photon reheating effect from entropy conservation, 
which was calculated earlier. Scaling this to the critical density Pc ~ ^Hq/SttG 
yields the relic's fraction of critical density. It is customary to define the 
measured value of the Hubble constant as Hq = 100 km s^^ Mpc^^. Then 



We can further reduce the expression for by noting that 

T^V s{T,) ^^Np{T^) 2 



T^J s{Td) Np{Td) Np 



(35) 



Using this relation and evaluating all numerical constants we get the convenient 
form 

^xh"" = , r^j, . (36) 
p,^ \/ N p J [x p) 

where /i = 1.2 x 10'"' GeV. The above formula is an accurate approximation to 
the Boltzmann equation solution to within about 15%. Table 1 lists the value 
of Np for various ranges of the freeze-out temperature Tp = xpm^. 

Since a massive weakly interacting particle decouples in the non-relativistic 
regime, one is able to expand the annihilation cross section in powers of the 
relative velocity 

av = a+^v'^ + ■■■ (37) 



The thermal average of this expansion yields 

{av) =a+ (b-^a\x + ---. (38) 
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With xp — 1/20 this is a quickly converging expansion. The constants a and 
h can be evaluated straight-forwardJ^ with knowledge of the squared-matrix 
element of the annihilation process iHl. Helicity amplitudes also exist for all 
MSSM processesB. 

The power series expansion of the thermal averaged cross-section breaks 
down, howeszer, when the relic particle annihilates into a pole (e.g., a XX ~^ ^ 
resonance) E^EJ. One then needs to use more careful techniques. Also, the 
Boltzmann equation becomes more complicated if there are other particles 
around with mass slightlx higher than the stable relic (within a "thermal dis- 
tance" |tox — < Tf)I1j. In this case, the close- by particles can help restore 
thermal equilibrium through coannihilation channels, effectively lowering the 
relic abundance. These potentially important cases will not be considered 
further here. 

4 Neutralino dark matter 

From the previous sections we have found that Q,^h? depends on l/{av). 
By dimensional analysis we can identify (av) ~ 1/™smsj/; which implies that 
oc rn^susy Since the age of the universe constraints arc compatible only 
with Q,^h? < 1 (very conservative requirement) then it should not surprise us 
that relic abundance considerations put a upper limit to how large rrisusy can 
be. We shall see quantitatively the results of these constraints in the following 
paragraphs. Before going straight to that, a few general comments about the 
lightest neutralino are useful to review. 

The neutralino is a majorana particle, meaning that the particle and con- 
jugate are the same. Therefore when two of them come together to-jannihi- 
late, Fermi statistics requires each to be in a different helicity state Ej. Thus 
XX ~^ f f requires a final state helicity flip for the external fermions, and the 
total annihilation rate in the s-wave is suppressed by to^/to^ compared to 
processes that do not require helicity flips. This is often referred to as "p-wave 
suppression" : the p-wave configuration doesn't have this Fermi statistics re- 
quirement, although it is suppressed by powers of w^. However, at higher values 
of rriy^ other channels start opening up, such as W^W~ and tt which do not 
have this suppression (although they may still have other coupling constant 
suppressions). 

The majorana nature of the neutralino is one important property of su- 
persymmetric dark matter. The other important realization is that numer- 
ous other supersymmetric particles play a role in the neutralino annihilation 
channels. The lightest neutralino, being the LSP, only annihilates into stan- 
dard model particles when the temperature is near freeze-out. However, many 
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important diagrams of these annihilation processes depend on intermediate 
supersymmetric states. For example, XX ~^ e~^e~ depends not only on an 
intermediate s-channel Z boson, but also on t-channel slepton exchange. Sim- 
ilarly, XX ~^ depends on t-channel squark exchange, and xX ~^ W~W~^ 
can depend crucially on the chargino spectrum. In general, the neutralino an- 
nihilation rate and therefore the LSP relic abundance depends on the entire 
supersymmetric low energy spectrum - this includes particle content, masses, 
and mixing angles. 

Upon closer inspection most choices of the supersymmetric spectrum ulti- 
mately give rise to an LSP relic abundance dependent on only a few parameters 
in the theory. Relic abundance in the minimal model depends most sensitively 
on only two parameters, the bino mass and the right-handed slepton mass. As 
discussed in Martin's introductory article in this volume, this minimal model is 
described partly by a common scalar mass mo at the high scale and a common 
gaugino mass mi/2- As long as mo is not much greater than toi/2 one typically 
finds in the low energy spectrum a bino (superpartner of the hypercharge gauge 
boson) as the lightest supersymmetric particle. One also typically finds the 
squarks significantly heavier than the sleptons, with the right-handed sleptons 
being lighter than the left-handed sleptons. 

Taking our hints from the minimal model, a useful initial exercise is to 
declare the LSP as a pure bino. Then, the largest contribution to the neutralino 
annihilation rate will be fmm t-channel Ir exchange in xX ^ l^l^ ■ This case 
has been studied in detail E3, where it was shown that 




(39) 



where M ~ 1 TeV and E — -I- m? . A good approximation to the require- 
ment that Q^h'^ < 1 yields 

mf 

— ^ < 200 GeV. (40) 
m^ 

This result has two important consequences. One, it satisfies the general obser- 
vation that the relic abundance constraint puts an upper limit on at least one 
superpartner. And second, in this important example of pure bino, both the 
bino and the right-handed lepton superpartner must have masses below about 
200 GeV. A priori the upper bound could have been any number. The result 
is compatible with weak scale supersymmetry {rrisusy ~ Tnw), long thought to 
be important in electroweak symmetry breaking. 
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It is not expected that the hghtest supersymmetric particle is a pure bino, 
but rather mostly bino. The coupling of the neutralino to the Z boson depends 
on the higgsino components. If the LSP is partly Higgsino then efficient an- 
nihilations through the Z boson could be possible, thereby reducing the relic 
abundance for a particular m-^ value. Limits on the superpartner spectrum still 
exist; they just happen to be at different scales than the 200 GeV we found for 
the pure bino case. 

Our example realistic particle spectrum is a precise formulation of the 
minimal model with electroweak symmetry breaking enforced. The full 4x4 
neutralino mass matrix is diagonalizejd numerically, and the relic abundance is 
calculated. The four contour plots of Fig. |l] demonstrate the effects of the 
relic abundance on the parameter space in the mo vs. mi/2 plane. The value 
of the top quark mass for this contour plot is 170 GeV and the sign of the 
supersymmetric /i parameter is chosen to be negative. In (a) tan/3 = 5, and 
^o/'TT-o — 0, (b) tan/3 = 5 and Aq/too — —2, in (c) tan/3 ~ 20 and A^jm^ — 3, 
and in (d) tan/3 = 10, Ao/nio = —2, and the sign of ^ has been changed to 
positive. The region inside the solid curve is allowed by all constraints. Letter 
labels are placed outside the solid curve to demonstrate which constraint makes 
the region phenomenologically unacceptable: A indicates Sl^/i^ > 1; B B{B — > 
sj) > 5.4 X 10"'*; C m + > 47 GeV; E electroweak symmetry breaking problems 
(full one loop effective potential becomes unbounded from below); and, L the 
LSP becomes charged (tij^ < rriy^). Current limits on rn^+ and B{b — > sj) 
are somewhat more restrictive than the limits quoted above, however they will 
only shrink the allowed parameter space a little bit near B and C. Inside the 
dotted line corresponds to one particular definition ofj^cceptable finetuning 
allowed in the electroweak symmetry breaking solutionsEEl. It is not important 
for our further discussion. 

The peak or spike in the allowed region for mi/2 — 130 GeV corresponds 
to XX annihilations through a Z pole. Since the LSP is not just a pure bino, 
and has some higgsino component to it, the Z pole annihilation rate becomes 
extremely important in this region. Also, it appears obvious from the figures 
that the relic abundance constraint has the most effect in the increasing niQ 
direction. This is understandable for several reasons. Away from the Z and 
Higgs poles, the LSP annihilations occur most efficiently through the t-channel 
scalar diagrams. For a fixed value, which is roughly equivalent to fixed 
TOi/2, the annihilation cross-section decreases rapidly as the scalar mass mo 
increases. Hence, a cutoff in the allowed region must occur as mo increases. 

Forays into the large mx/2 region with mo fixed usually end up with other 
problems, the most common of which is Iji becoming the LSP. This is sim- 
ply a result of the renormalization group equation dependences on the U{1)y 
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Figure 1: Contours of allowed parameter space, mt = 170 GeV for each frame. In frame 
(a) tan/3 = 5, sgn(/i) = — , and Ao/mo = 0, (b) tanfi = 5, sgn(/i) = — and Ao/mo = —2, 
in (c) tan/3 = 20, sgn(/n) = — and Ao/mo = 3, and in (d) tan ,8 = 10, sgn(/i) = +, and 
Ao/mo = —2. See the text for explanation of letter labels. 
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couplings. A quick inspection of the renormalization group equations for to^ 
and for mj^ proves that if mi/2 3> mo,mw, then the right-handed slepton 
must be hghter than the bino. RuHng out this region of parameter space can 
be considered a_reHc abundance constraint, since the dark matter probably 



Realization that a weak-scale higgsino could be a legitimate dark matter 
particle is a rather recent development. One way to obtain an higgsino as the 
lightest neutralino is to make much less than the gaugino parameters in the 
neutralino mass matrix. A very low value of /x will create a roughly degenerate 
triplet of higgsinos. The charged higgsino and the neutral higgsinos can all 
coannihilate together with full SU{2) strength, allowing the LSP to stay in 
thermal contact with the photons more effectively, thereby lowering the relic 
abundance of the higgsino LSP to an insignificant level. These coannihilation 
channels are often cited as the reason why higgsinos are not viable dark matter 
candidates. This claim is true in general, but there are two specific cases that 
I would like to summarize below that allow the higgsino to be a good dark 
matter candidate. 

Drees et al. have pointed out that potentially large oneJjpop splittings 
among the higgsinos can render the coannihilations less relevantliS. Under some 
conditions with light top squark masses, one-loop corrections to the neutralino 
mass matrix will split the otherwise degenerate higgsinos. If the mass difference 
can be more than about 5% of the LSP mass, then the LSP will decouple from 
the photons alone and not with its other higgsino partners, thereby increasing 
its relic abundance. 

Another possibility lil relating to a higgsino LSP is to set equal the bino 
and wino mass to approximately mz- Then set the — /i term to less than mw- 
This non-universality among the gauginos and particular choice for the higgsino 
mass parameter, produces a light higgsino with mass approximately equal to 
fi, a photino with mass at about mz, and the rest of the neutralinos and 
both charginos with mass above myy- There are no coannihilation channels to 
worry about with this higgsino dark matter candidate since no other chargino 
or neutralino mass is near it. The value of tan/? is also required to be near 
one so that the lightest neutralino is an almost pure symmetric combination 
of Hu and iJ^ higgsino states. The exactly symmetric combination does not 
couple to Z boson (at tree level) . The annihilation cross section near tan /3 ~ 1 
is proportional to cos^ 2/3. The relic abundance scales inversely proportional 
to this, and so the nearly symmetric higgsino in this case is a very good dark 
matter candidate. Note that there are no i-channel slepton or squark diagrams 
since higgsinos couple to sfermions proportional to the fermion mass. Because 
the higgsino mass is below mw, the top quark final state is kinematically 
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inaccessible, and so the large top Yukawa cannot play a direct role in the 
higgsino annihilations. 

This non-minimal higgsino dark matter candidate described in the pre- 
vious paragraph was motpated by the e"'"e~77 event reported by the CDF 
collaboration at FermilabEd. The non- minimal parameters oil which leads to 
a radiative decay of the second lightest neutralino (photino) into the lightest 
neutralino (symmetric higgsino) and photon also miraculously yield a model 
with a good higgsino dark matter candidate. 



5 Conclusion 

The minimal model bino and the higgsino described above work as dark matter 
candidates both qualitatively and quantitatively. Nature, of course, might not 
conform to either of these specific possibilities, but it is straight-forward to 



'here are, of course, numerous ways 
Scalar mass universality could be 



catalog the non-minimal possibilities. 

to go beyond what is presented here I 

relaxedo. One could in fact suppose that i?-parity is not exactly conserved 
and the LSP is not a stable particle. Perhaps one could allow small i?-parity 
violations which create meta-stable LSPs whose lifetimes are greater than the 
age of the universe to solve the dark matter dilemma. However, it is not enough 
to just make the lifetime greater than the age of the universe. Remarkably, the 
measurements of positrons and photons in cosmic rays require that the LSP 
lifetimes into-these particles be many orders of magnitude beyond the age of 
the universe o. In other words, it is not easy to make the LSP a meta-stable 
dark matter candidate. If i?-parity were broken in nature, it appears necessary 
to look elsewhere for the dark matter candidate. 

One common theme exists in non-minimal models which accommodate a 
supersymmetric solution to the dark matter problem. It is the requirement that 
XX not couple to the Z boson. A full strength SU{2) coupling to the Z boson 
generally allows too efficient LSP annihilation, and therefore too small relic 
abundance to be relevant to the dark matter problem. This theme ran be found 
in several noiL^ininimal examples such as the symmetric HiggsinoE£IE3 described 
above, the Zl3, and the sterile neutralinoc3 dark matter candidates. Of course, 
the minimal model also provides a low-strength coupling to the Z naturally 
with the mostly bino dark matter candidate lZI. The older supersymmetry dark 
matter literature considered the photinii, which also does not couple to the Z, 
as the primary dark matter candidate E3. 

Theoretical ideas about the low-energy superpartner spectrum always have 
dark matter implications. The supersymmetric solution to the dark matter 
problem is perhaps only superseded by gauge coupling unification in exper- 



15 



imcntally based indications that nature might be described by softly broken 
supersymmetry. It is mainly for this reason that pauch experimental effort is 
expended on the search for supersymmetric relics c3. Likewise, theoretical ef- 
forts on understanding the origin of LSP stability and the prediction of dark 
matter properties should continue to enlighten us. 
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